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We show that rather simple but non-trivial boundary conditions could induce the appearance of spatial chaos 
(that is stationary, stable, but spatially disordered configurations) in extended dynamical systems with very sim- 
ple dynamics. We exemplify the phenomenon with a nonlinear reaction-diffusion equation in a two-dimensional 
undulated domain. Concepts from the theory of dynamical systems, and a transverse-single-mode approxima- 
tion are used to describe the spatially chaotic structures. 
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I. INTRODUCTION 



In the past few decades, considerable understanding of the 
phenomenon of temporal chaos in dynamical systems of few 
degrees of freedom has been achieved [|l]-p|]. On the other 
hand, spatio-temporal chaos in extended dynamical systems 
with infinitely many degrees of freedom is currently under 
very active investigation [|[^]. It is remarkable however, that 
an area of problems laying somehow between the two ex- 
tremes has not received so much attention, namely, purely 
spatial chaos as a stationary attractor of extended dynamical 
systems |||-[l^]. 

The possible existence of this kind of attractors was first 
suggested by Ruelle [ 1 1 1 in the context of equilibrium phases. 
He pointed out the parallelism between a time-dependent dif- 
ferentiable dynamical system and the space dependence of 
equilibrium states in statistical mechanics. He then raised 
the question as to whether the existence of turbulent crys- 
tals could be the natural next step towards complexity after 
spatially homogeneous, periodic, and quasi-periodic equilib- 
rium phases have been found. Newell and Pomeau [|l2|| gave 
some conditions under which such a turbulent crystal would 
exist in pattern-forming systems described by a free energy. 
Theoretical and experimental work on modulated phases and 
commensurate-incommensurate transitions [[Ij] represent ad- 
ditional concrete results along these lines. 

In the context of fluid dynamics, the existence of spatially 
chaotic, but temporally steady solutions would also fill a con- 
ceptual gap between two well-studied complex phenomena: 
Lagrangian chaos, and Eulerian chaos or turbulence. The 
former refers to the chaotic motion of a fluid parcel which 
might occur even in laminar and, in three dimensions, steady 
flows Jl4|4l6|]. On the other extreme, the road to turbulence is 
usually associated to a hierarchy of instabilities leading to in- 
creasingly spatio-temporally chaotic Eulerian velocity fields. 
Frozen spatial chaos would then refer in this context to a third 
possibility: a stationary flow spatially chaotic in the Euler de- 
scription. 

By now, many extended dynamical systems displaying spa- 
tial chaos have been identified. Most of the studies are con- 
cerned with one-dimensionally extended systems. They are 
specially suitable to analysis because their steady state config- 



urations depend just on the unique spatial coordinate. These 
configurations are solutions of sets of ordinary differential 
equations (the spatial dynamical system) with the space vari- 
able as the independent variable. The standard theory of low- 
dimensional dynamical systems can be used to describe such 
configurations, by just considering the spatial coordinate as a 
fictitious time. Rigidly traveling waves with spatial chaotic 
structure can also be considered as a case of spatial chaos in a 
moving reference frame ||^,^ . 

Spatial chaos may appear when the spatial dynamical sys- 
tem has a sufficiently high dimensional phase space. This 
high dimensionality may arise from either a) the presence of 
high-order spatial derivatives in a single evolution equation 
as in the cases of the Swift-Hohenberg equation 1 17,181], and 
Kuramoto-Sivashinsky and related models [pHlO|], b) the cou- 
pling of several fields each one satisfying a lower order differ- 
ential equation as in (the real and imaginary parts of) the com- 
plex Ginzburg-Landau equation [|| which supports chaotic 
travelling waves, or c) explicit space dependent forcing terms 
as in 1] or ^ . 

Consideration of two-dimensional spatial chaos has been 
much scarce. In absence of a simple connection with conven- 
tional dynamical systems theory, the very concept of chaos 
in the spatial configuration should be properly defined for 
the general two-dimensional case. A rather complete formal- 
ism generalizing dynamical system tools (entropies, dimen- 
sions,.. J to multidimensional spatial chaos has been devel- 
oped [[7pC|], and some examples examined [pljlq]. On the 
other hand, as the number of dimensions increases, a much 
larger variety of non-trivial boundary condition classes surely 
leads to a greater richness in the expected properties of the 
steady field configurations. A well-posed question is then 
whether relatively simple boundary conditions may lead to 
steady spatially chaotic configurations. The main purpose of 
this Paper is to address this question. 

In addition to the existence of chaotic spatial configura- 
tions, it is important to study also their stability in time. A 
stationary state will only be physically observable if it is stable 
or at least long-lived. It turns out that the temporal stability of 
the stationary solutions is in general unrelated to the stability 
of these configurations considered as orbits of the spatial dy- 
namical system. In the examples cited above, there are cases 



of both stable and unstable space-chaotic configurations, but 
instability seems to be more frequent. As a consequence, spa- 
tial chaos has been generally considered of limited physical 
relevance. 

In this Paper we show that rather simple undulated strip- 
like domain shapes can induce, in a very simple nonlinear ex- 
tended dynamical system, the formation of patterns that are 
both spatially chaotic and temporally attracting. The kind of 
modulated boundaries we use could be easily implemented 
in standard experimental pattern-formation set-up's such as 
Faraday waves, convection cells, or open flows. In fact, our 
work was originally motivated by the observation, in a fluid 
dynamics experimental setup consisting of a periodic array of 
pipe bends, that the transverse profile of the steady flow does 
not necessarily repeat itself with the same periodicity of the 
array ^. 

In Section y we present the particular model that we study 
and perform a preliminary analysis of its behavior. In Sec- 
tion |l| a single-transverse-mode approximation is introduced 
and we use it to predict the existence of boundary-induced 
spatial chaos. Numerical simulations are presented in Sec- 
tion 1^ to substantiate our claims beyond the validity of the 
previous approximation. Finally we summarize the results and 
open problems in the Conclusions. 



II. A REACTION-DIFFUSION EQUATION IN A 
STRIP-SHAPED DOMAIN 

As stated before, the application of the theory of dynam- 
ical systems to the study of stationary spatial configurations 
of one-dimensionally extended systems is direct. A stationary 
pattern satisfies, in general, a system of ordinary differential 
equations with the spatial coordinate as its independent vari- 
able which we can think of as a time. Parity symmetry in a 
spatial coordinate will appear as time-reversal symmetry after 
reinterpretation of this coordinate as time. 

The general study of spatial chaos in several spatial dimen- 
sions requires the notion of translational dynamical systems 
with d times [0,|o|]. There are situations, however, where such 
formalism is not necessary because the two independent spa- 
tial directions are distinguished by the geometry of the system, 
so that one of them naturally plays the role of time. In this 
way, the spatial variation in one direction would be interpreted 
as time evolution of a one-dimensional field that only depends 
on the remaining spatial coordinate. Particularly suited to our 
approach will be the case of two-dimensional extended sys- 
tems in strip-shaped regions much longer (ideally infinite) in 
the time-like direction than in the space-like one. If the strip is 
narrow enough, only patterns composed of one or few trans- 
verse spatial modes will be allowed and spatial chaos could be 
readily defined and identified in terms of the usual concepts of 
dynamical systems theory. 

In order to concentrate on spatial chaos purely induced 
by boundary effects, we consider a very simple model equa- 
tion containing only up to second order derivatives and a sin- 
gle field variable, a reaction-diffusion equation of the Fisher- 
Kolmogorov type: 



dtip = V'^tp + aip-ip^ 



(2.1) 



with appropriate boundary conditions for the real field ip{x, t). 
The real linear coefficient a can be absorbed rescaling the 
variables, but we find convenient to keep it explicit in the 
equation. Equation {TA} appears in several contexts including 
phase transitions, where it takes the name of real Ginzburg- 
Landau equation or time-dependent Ginzburg-Landau model 
pi]], and population dynamics [p2|]. The dynamics of (2.1) 
can be written as purely relaxational [ P3| ] in a functional Lya- 
punov potential V[i/j]: 



with 
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D 



dt^ 



dx ( -IWp 
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(2.2) 



r + 7^ + sm ' (2-3) 



where the integral is over the domain D. The surface term 
Sl'ip] takes into account the effects of boundary conditions 
over the domain limit, and it vanishes when periodic, null 
Dirichlet (-0 = 0) or null Neumann {dni} = 0) boundary con- 
ditions are specified. It follows from Eq. ( |2.2| ) that V can only 
decrease with time. The relaxational character of Eq. ( |2.2| ) 
implies also that the only asymptotic states are fixed points. 
Therefore, this model does not display any limit-cycle oscil- 
lations or more complex dynamics such as temporal chaos in 
any number of sp atial dimensions. 

Equation (2.1) has been extensively studied in one and two 
dimensions. In one dimension, for infinite systems, we have 
the following situations: For a < 0, i/" = is the only sta- 
tionary solution and is stable under time evolution. At a = 
a pitchfork bifurcation occurs and the former trivial solution 
looses its stability. For a > some of the stationary solutions 
are the following: 

1. Homogeneous solutions: ip{x,t) — 0, and i(;{x,t) — 
±^/a = tI)±. 

2. Kink type solutions: %lj{x^t) = ±-/atanh(w^a;). 



3. Periodic solutions: ■tp{x,t) = y/a — k^ sin(fca;) + ..., 
where the ellipsis stands for higher order harmonics. 

Only ■tp± and the kinks are now linearly stable. 

For finite systems, the boundary conditions normally re- 
strict the variety of solutions either by selecting some values 
of k, or by favoring either the kink or the ijj± solutions. Ad- 
ditionally, the boundary conditions may have as an effect a 
shift in the pitchfork bifurcation point a = to a different 
value flc 7^ and a change in the coefficients appearing in the 
solutions. 

The typical time evolution of initial field distributions leads 
to the formation of domains where values of ip close to either 
ijj+ or V'- dominate. These domains are separated by kink- or 
anti-kink-type walls that can move into each other producing 
a mutual annihilation. By this mechanism, small domains dis- 
appear and feed the larger domains whose sizes then increase 



logarithmically in time until one of the stationary solutions, 
prevailing by chance, takes over the whole system. 

In two dimensions, the dynamics typically consists of the 
coarsening of domains of if)^ and ip- phases whose typical 
size grows as the square root of time |21 1. Interfaces between 



the two phases are locally similar to t he on e-dimensional kink 
solutions. The gradient term in Eq. (2.3) is important at the 
interfaces giving a positive contribution to the Lyapunov po- 
tential. Since the dynamics always minimizes V, it tends 
to reduce the length of these interfaces. This reduction is 
achieved by the shrinking and ulterior collapse of the small- 
est domains to contribute to the (square-root) growth of the 
remaining ones. 

It was shown by Collet p4| ] in a more general context that 
the time evolution and final states of Eq. ( 2J^ ) in finite domains 
are similar to those in an infinite system except in a boundary 
layer around the border whose size depends on the a parame- 
ter. This result holds in both one and two dimensions. Thus, 
in order to observe the influence of boundaries on pattern evo- 
lution, we need to consider a domain small enough at least 
in one of the directions. In a stripped domain, elongated in 
the X direction, this small dimension will be the transverse y 
direction. The domain will be hmited in this transverse di- 
mension by the boundaries yQ{x) and yi{x), where the func- 
tion i/j{x,y,t) will take values iIjq{x) y ipi{x), respectively 
(Dirichlet conditions). 

In the one-dimensional case, conventional dynamical sys- 
tems theory implies that there are no chaotic stationary solu- 
tions to Eq. (Z.l) because the spatial dynamical system is just 



a second order ordinary differential equation. However, chaos 
can arise if some a;-dependent periodic forcing is added to the 
equation. These arguments do not apply directly to the two- 
dimensional case. However, it is tempting to think of undu- 
lations of the lateral boundaries as a kind of periodic forcing 
on the longitudinal coordinate. This suggests the possibility 
of finding chaotic structures in the x-direction induced by un- 
dulated boundaries. 





FIG. 1. Examples of stripped channels enclosed in oscillating 
walls: a = 27r/A; do = rfi, and (a) (p — n , (b) tf) — 0. 



As a particular case we consider domains limited by two 
sinusoidal boundaries. With applications to hydrodynamics in 
mind we think of these domains as channels with sinusoidal 
banks. 



yii^) ^ y(1 ~ cos(aa;)) (2.4) 

yo(a;) = -l-y(l + cos(aa; + 0)) (2.5) 

Here di, do are the amplitude of the undulation of each 
bank, a is the spatial frequency which we assume to be the 
same for both banks and (p is their mutual phase mismatch. 
Figure (p shows a few typical shapes for our channel-like do- 
mains. The case di — do and (f) — n gives a sausage shaped 
channel with symmetrically and sinusoidally varying width. 
On the other hand, di = do and = sets the boundaries in 
phase and corresponds to a domain with the form of a sinu- 
soidally meandering channel of constant y-width. 

We stress that we want to consider the simplest situa- 
tion that may display spatial chaos. Consideration of more 
complex equations exhibiting spatial chaos even with simple 
boundaries, or more complex boundaries such as incommen- 
surate oscillations for the upper and lower banks (correspond- 
ing to quasi-periodic forcing) would only enrich the complex- 
ity of stationary solutions. 

A complete definition of the model requires also the specifi- 
cation of boundary conditions on the longitudinal x direction. 
The analogy with a temporal variable would be better for do- 
mains infinite in the x direction, with only the weak require- 
ment of boundedness for ip. However, an infinite domain is 
inadequate for the numerical approaches to be described be- 
low. In our calculations we would need to impose periodic 
boundary conditions (of period L) along the x direction. In 
this way we are restricting the class of solutions to periodic 
orbits of period L or less in the time-like coordinate. We 
will still be able to identify as spatially chaotic the configu- 
rations that have the maximal period L, provided this period 
increases and the periodic orbit approaches a chaotic trajec- 
tory as system size L increases. Subtle considerations such as 
Lyapunov number computations for such limiting orbits will 
be addressed elsewhere. 

To perform numerical simulations we can choose between 
several strategies. If we are only interested in stationary states, 
we can numerically solve the time-independent version of 
Eq. (2.1) by means of finite elements or finite differences. 



These methods can be implemented to find solutions that may 
or may not be stable under time evolution. Another possibility 
is to follow the dynamics of the full Eq. (2.1 ) until a station- 
ary state is reached. In this way, only attracting (i.e. stable) 
stationary solutions can be found (remind that only stationary 
attractors are allowed by this purely relaxational dynamics). 
In any case, a convenient way to handle the boundary con- 
ditions is to map the region limited by yo (x) and j/i (x) (and 
by x = 0, L) to a rectangular one: j/i — 1, yo = 0, and 
X ^ 0,L. For arbitrary functions yo{x) and yi{x), the map 
{x,y) I — > (x,2/),with 



y^iy- yo)/{yi - yo) , 



(2.6) 



transforms Eq. (2.1) into an equation for ■>p{x,y,t) 
ip{x,y,ty. 



where 






g{x) = -2 



A2 
^xV + Vox 



nix 

A{x 
Ax{x 



y{2Al - AA:,:,) - Ayoxx + 2A^yox 



M^) = yi{x) -yo{x) 

A{x) 



dx 

(f 
Axx{x) = -7-^A{x) 

d 
yo,ix = 3-2/0,1 

dx 



(2.7) 

(2.8) 

(2.9) 

(2.10) 
(2.11) 
(2.12) 

(2.13) 
(2.14) 



Here A(x) is the transverse distance between the boundaries 
and gives the width of the strip. If this width does not vary 
along x, then Ax{x) = = A2;a;(x). 

The new transverse boundary conditions are 



^{x,y = 0)='ijjo{x), ipix^y = 1) = ipi{x) 



(2.15) 



A first observation is that the shape of the domain bound- 
aries is reflected as a parametric forcing of the equation in the 
new coordinates. For example, fo r th e simplest case of a me- 
andering channel given by Eqs. (2.4) and (2.5) with = 0, 
and do = di = d (Fig. [lib): 

dti' = dlJ + H^)dy: + g{x)dl^^ 

+ nix)dyrp + aiP-iP^ (2.16) 



Tix)^ 
g{x)^ 
H{x) = 



da sin{ax) 
2 



(l + d)2 

da sin(Q;a;) 

l + d 
da^ cos(aa:) 



2(1 + d) 



(2.17) 
(2.18) 
(2.19) 



Setting the right hand side of Eq. (2.7) to zero in order to 
seek for statio nary solutions, and thinking of x as the time we 
can view Eq. ( Z7 ) as a nonlinear evolution equation for a one- 
dimensional field with a 'time' -periodic parametric driving 
due to the boundaries. Present knowledge on spatio-temporal 
chaos and pattern formation can in principle be applied to ana- 
lyze the behavior of this resulting evolution equation. General 
results are not abundant, however In the next Section fur- 
ther approximations will be introduced in order to facilitate 
the analysis and establish the existence of stationary spatially- 
chaotic solutions. 



III. SINGLE-TRANSVERSE-MODE APPROXIMATION 

For definiteness, in the rest of the Paper we will consider 
just null Dirichlet boundary conditions, that is the field ip 
takes the value zero at the transverse boundaries: tpojx ) — 
i/jiix) = 0. In this case our model in the form of Eq. (2.7) 
has the trivial solution ilj{x,y,t) = 0. The stability anal- 
ysis of this solution for the case of a rectangular domain 
of width I leads to an eigenvalue problem for the linearized 
equation. For a > Uc — (j)^ the eigenfunctions factor- 
ize into longitudinal (exp (ikxx)) and transverse (sm{kyy)) 
modes: ip\,k{x,y,t) = exp{Xt) exp{ikxx)sm{kyy), with 
ky = jm,m — 1,2,..., k^ real, A satisfying the disper- 
sion relation A = a ~ k^, and k"^ = kl + kj. The unsta- 

X y 

ble modes are then those satisfying the condition k^ < a. 
The first unstable mode corresponds to {kx,ky) — (0, j), 
which becomes unstable at the critical value Qc — ( j)^- The 
transverse modes are discretized in multiples of j due to the 
boundary condition. If in addition, we require L-periodicity 
in the longitudinal coordinate x, k^ will be also discrete, but 
provided that L ^ / this discretization will be much finer 
than the transverse one. The value of the parameter a con- 
trols how many modes are linearly unstable. If the trans- 
verse size I is small enough for the control parameter to sat- 
isfy the condition (j)^ < a < (^)^, there would be just 
one linearly unstable transverse mode, with many associated 
longitudinal unstable modes. Close enough to the instabil- 
ity threshold, we can try an approximate solution of the form 
i^ix, y, t) = A{x, t) sin(7ry) and write an evolution equation 
for the amplitude A(x, t) of the first transverse mode. 

We are interested however in a domain which is not a rect- 
angular strip but a undulated channel. The coordinate change 
bringing our domain into a rectangular one renders the vari- 
ables non-separable and the linear problem is no longer solv- 
able analytically. However, for small deviations from the uni- 
form channel a perturbation scheme can be used. In the same 
vein as in the preceding paragraph we try, for our undulated 
domain, the ansatz ■0(2:, y, t) = A(x, t) sin(7ry) + O{do, di), 
assuming that a is close to the threshold imposed by the 
small I, and that the size of the channel undulations is small. 
We call this approach a single-transverse-mode approximation 
(STMA). 



Projecting Eq. (2.7 ) onto the single transverse mode present 
in the ansatz and neglecting higher order contributions we get 
the following evolution equation for the amplitude A: 



dtA = dl^A + Pix)dxA + uj^{x)A - -A^ (3.1) 



where 



Pix) = -^ (3.2) 



and 



2, - AA^^-2A| /7r\2 



(27r2 



(3.3) 



(a) 



We have checked that a more rigorous, but lengthy, ap- 
proach based on a multiple scale expansion leads to the same 
result [pSl. The stationary patterns satisfy the time indepen- 
dent version of Eq. (3.1). In terms of the coordinate x con- 
sidered as a time, one can view this spatial dynamical system 
as a parametrically forced nonlinear oscillator: domain un- 
dulations provide a periodic driving on the frequency of the 
system tj^(a;). In turn, /3(x) is a "dissipation" term which can 
be positive or negative depending on x. This x-modulation 
of (3 comes from the longitudinal variation of the vertical 
width. The integral of the "dissipation" (3{x) on one period 

of the oscillating boundaries T = 2Tr/a is / " (3{x)dx = 



In 



A(x+T) ~ ^- "^^"^^ shows that although system ( |3.1[ ) is 
locally dissipative, it is effectively conservative over one pe- 
riod of the modulation. This implies, in particular, that the 
stroboscopic map associated with the system (3.1 ) is area pre- 
serving. 



In the general case, Eq. (3.1) can be simplified by re- 
moving the dissipation term with the change: A{x, t) = 
exp (— i J/3(a;)(ix) p(x,t) = p{x , t) / ^/A{x) . The new 
equation reads: 



dtp^dl,p + n\x)p~—p' 

(-TT \ 2 
^j (1 + 2/ix2/ox) - 



12 



-0' 



(3.4) 
(3.5) 



A particular case occurs when the transverse distance be- 
tween the two channel borders does not vary along the longi- 
tudinal direction, so that the dissipation term vanishes identi- 
cally, i. e. P{x) — 0. For the sinusoidal channel this happens 
when = and do = di = d (Fig. [l|b). In this case, the 
amplitude equation is reduced to 



dtA = dl,A + u;\x}A~-A^ 



(3.6) 



with 



(nda) 
4(1 + d): 



■ sin (ax) 



(3.7) 



When dfA ~ 0, (3.6) is a dissipative Mathieu equation modi- 
fied by the addition of a cubic nonlinear term or, equivalently, 
a parametrically forced Duffing oscillator This equation is 
known to have chaotic solutions \UM . 
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FIG. 2; ^(a) Stroboscopic Poincare map of the phase space of the 

system (3.8) for the values do = ^i = 1. ^ = 27r, a = 17, (f) = n. 
KAM tori and chaotic trajectories inbetween are clearly seen. The 
arrow indicates the approximate location of the fractal boundary sep- 
arating bounded and unbounded trajectories, (b) The chaotic config- 
uration corresponding to the cloud of points surrounding the origin 
in the stroboscopic map. 



The general time-independent case of Eq. (3.1) [or 
Eq. K^] reads: 



d^^A = -iu^{x)A + -A^ - f3{x)d.,A . 



(3.8) 



In the absence of undulations (do — di = 0) paramet- 
ric forcing and dissipation vanish and the equation is both 
Hamiltonian and integrable. It has a region in phase space 
{A, Ax) close to the origin where motion is bounded and 
regular Beyond the separatrices of the two saddle points 
{A, Ax) = ±2/\/3(a — 7r^,0), trajectories escape to infin- 
ity. When undulations are introduced, separatrices of the sad- 
dle points deform and may cross. It is well known that for 
perturbed Hamiltonian systems, separatrix intersections indi- 
cate the onset of chaos. In our system, in addition to chaotic 
bounded trajectories, separatrix intersections lead also to frac- 
talization of the phase-space boundary dividing bounded and 
unbounded trajectories. Melnikov theory provide us with the 
tools to determine analytically the necessary conditions for 



separatrix intersection and the occurrence of chaos. Follow- 
ing [26 1, the Melnikov function M{9) can be calculated for 
small d. For example, in the case (/) — tt, do = di = d one 
finds: 



M{0)^ f{a,a)smiae) 



(3.9) 



The fact that the function M{d) has zeros as a function of 
6 indicates that separatrix intersection and chaotic behavior 
occur. Chaotic behavior appears for arbitrarily small sizes of 
the channel undulations. We expect similar behavior to occur 
for other values of ip, do and di . 
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FIG. 3. Stationary solution of the STMA ( |3.l| ) compared to a 
fully two-dimensional simulation of (2.1) for the parameter values 
do = rfi = 0.1, a = TT, a = 10, </> = n. a) Stationary solu- 
tion A{x) of the STMA (dotted) and the on-axis values of the actual 
two-dimensional stationary solution ip{x,y = 0) (solid), b) The 
full two-dimensional solution of (2.1) represented on a gray-scale. 



White corresponds to the highest values of ip and black to the low- 
est, c) Reconstruction of the two-dimensional field from the STMA 
{tp(x, y) — A{x) sin(7r?y)) on the same gray-scale as in b). 

To illustrate the chaotic behavior of the stationary STMA 
(3^) we show in Fig. (ga) its numerical stroboscopic map. 
The dots are values of {A, A^) at multiples of the forcing pe- 
riod T — 2TT/a for a set of initial conditions. Several regions 
dominated by chaotic trajectories and separated by KAM tori 
(the closed curves corresponding to quasi-periodic solutions) 
are clearly recognized in the picture. Also, the approximate 
location of the fractal boundary separating bounded trajecto- 
ries from those escaping to infinity is pointed by an arrow. 
Melnikov analysis also implies the existence of a dense set of 



unstable periodic orbits in the vicinity of the separatrix inter- 
section on both sides of the fractal boundary mentioned be- 
fore. These periodic orbits and the fact that they constitute 
a skeleton of the chaotic trajectories, are important for our 
analysis because the periodic boundary conditions in the a;- 
direction select them out from the uncountable many other 
possible solutions of Eq. (3.8). 



IV. NUMERICAL STMA AND TWO-DIMENSIONAL 
TIME-INTEGRATIONS 



In the previous section, we have shown both analitically 
and numerically that simple undulated boundaries may induce 
spatially chaotic steady solutions in our simple model (3.1). 
In this section, we discuss the accuracy and range of validity 
of the STMA and the physical relevance of its solutions by 
comparing with the numerical integration of the full model 
(P). 



In the first place, we discuss the results of such integra- 
tion for the particular case (p = tt and d^ = di = d. Start- 
ing from random initial conditions the system, after a time 
long enough, settles in a disordered stationary configuration, 
as shown in Fig. (g|). There we present the stationary con- 
figuration A(x) = A{x,t — > cx)) obtained by direct integra- 
tion of Eq. ( |3.l[ ). In the same plot, a longitudinal cross sec- 
tion of the asymptotic field obtained from a simulation with 
the full two-dimensional model is shown for comparison. We 
use as the initial condition of the two-dimensional problem 
the solution from the STMA for the same parameter values, 
■0 = A{x) sin(7ry), to provide an approximate stationary so- 
lution to Eq. (2.1). After a short time of adjustment, the sys- 
tem settles in a stationary state that is very close to the initial 
approximation. The full two-dimensional field and its recon- 
struction from the STMA are also shown in Figs, (gp) and 
(He), respectively. This figure reveals a strikingly accurate 
fit of the STMA solution and the complete field simulation, 
a strong indication of the validity of the approximation as a 
tool for analysis. The maximum absolute error of the approx- 
imate solution is of the same order of magnitude in both the 
undulated channel and a rectangular-domain test-case. This 
suggests that the error is mainly due to the truncation at the 
first linear transverse mode, but not to the peculiarities of the 
curved boundaries. 

The accuracy of the STMA, A{x, t) suviiry), breaks down 
when the channel width increases or when strongly non rect- 
angular domains are considered. Nevertheless, we have also 



performed direct simulations of Eq. (2.1) for this last case. An 
example of the typical behavior is shown Fig. (Q) for bound- 
aries defined hy d^ — di = 1.0, a = tt, a = 20, ip — tt. 
Notice that the resulting stationary configuration displays the 
same qualitative features of the STMA solutions: disordered 
distribution of kinks randomly pinned at some of the narrows 
of the channel. This is an evidence of the fact that undulating 
boundaries may also be the source of stationary spatial chaos 
in a two-dimensional system (2.1 ) beyond the regime well de- 



scribed by a single transversal mode. 
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FIG. 4. Two-dimensional steady state obtained by simulation of 
the system (2.1) starting from random initial conditions. The ampli- 



tude of the field on the channel axis is shown in the lowest panel. 
Parameter values: d = 1.0, a = tt, a = 20, (p = ir. 



The results in Figs. (0) and (0) illustrate the physical mech- 
anism behind the emergen ce o f spatial chaos in our system. 
Let us remember that Eq. ( |2.1[ ) evolves to minimize the po- 
tential (2.3). This minimization requires to reduce as much as 
possible the length of the interfaces between jp+ and tp- . Fol- 
lowing this tendency, an interface that links the opossite lateral 
banks of the channel and is far from any other interface will 
evolve to lock into one of the narrows of the channel, where it 
is shorter than in any other position. Detaching the interface 
from the bank of the channel would imply a temporary in- 
crease of the potential V due to the necesary proliferation of 
new interfaces. Such a potential increase is not allowed by the 
dynamics. Hence, the random occupation (arising from ran- 
dom initial conditions) of the narrows of the channel by kinks 
and anti-kinks finally builds up a spatially chaotic stationary 
configuration. This argument, based on kink and interface 
dynamics, clearly applies beyond the range of validity of the 
STMA, where we have analytically established the existence 
of spatial chaos, as evidenced in Fig. (Q). With this mechanism 
in mind, we can also conclude that not all the chaotic trajec- 
tories presented in Fig. (0) will lead to spatial chaos stable in 
time: only those corresponding to an energetically favorable 
(a local minima of the potential V) distribution of kinks will 
be reached under time evolution. In particular the trajectory 
plotted in Fig. (0b) corresponds to a temporally unstable con- 
figuration. 
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FIG. 5. One-dimensional solution of the system ( p.l\ ) showing a 
steady chaotic configuration, its spatial power spectrum, the trajec- 
tory in the projected phase space {A, Ax), and the time evolution of 
the Lyapunov functional, showing that this is an attracting configu- 
ration. Parameter values as in Fig. (Bl), but for a longer system 



Let us now try to get further insight about the chaotic nature 
of the irregular spatial structures described above. Periodic 
boundary conditions in the longitudinal direction always force 
the system to converge not to a chaotic spatial configuration, 
but to a periodic one wich we have shown it may very well 
be of the maximal period. To justify the use of the chaotic 
qualifier we need to show that as the size of the system in- 
creases, these periodic configurations approach one that could 
be characterized as chaotic in some way. Of course, to numer- 
ically carry out this process we would need to consider very 
long channels. Unfortunately though, performing direct sim- 
ulations on the fully two-dimensional model soon becomes 
computationally prohibitive as the number of the channel un- 
dulations increases. However, having demonstrated that the 
STMA accurately describes the qualitative features of the full 
model, we can concentrate our attention on the behavior of the 



approximate model Eq. (3.1 ) 



In Fig. (Q) we summarize the results from the numerical 



integration of Eq. (3.1) using the same parameter values as 
in Fig. (|^) with the exception of the domain size which now 
is much larger. An asymptotically stable configuration is 
shown in a) while c) displays the projection of the trajectory 
in phase space {A, A^)- The power spectrum of A{x) plotted 
in Fig. (pp) shows the typical broadband feature characteristic 
of chaotic trajectories. As a further indication of the (aproxi- 
mately) chaotic nature of this configuration, we have shown in 
Fig. (q) the stroboscopic map constructed from the trajectory 
in Fig. (ph. taking phase space points at 'times' integer mul- 
tiples of T. This map reveals an incipient self-similar frac- 
tal structure, also common in chaotic trajectories. All these 
facts together give compealing evidence that the trajectory, al- 
though L-periodic by construction, develops chaotic features 
as system size increases. 
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FIG. 6. Stroboscopic map of the trajectory shown in Fig. {|5|). 
Two succesive amplifications are shown, displaying a self-similar 
structure typical of a fractal. 



terns commonly used to study pattern formation, provided 
that boundary conditions such as those studied here are im- 
posed. The stability of these chaotic configurations should 
be discussed in each particular case. While in our model 
stabili ty c omes from the tendency of the dynamics to mini- 
mize ( Z3 ) therefore minimizing interface lengths and leading 
to pinning of these interfaces to the narrows of the channel, 
the mechanism for stability in other systems may be differ- 
ent. Apart from the direct application to pattern forming sys- 
tems, the idea of boundary-generated spatial chaos could be 
speculatively transferred to other nonlinear extended dynam- 
ical systems of interest. For example, it is possible that low 
Reynolds number fluid flows through a space-periodically per- 
turbed pipeline or even through a realistic channel of shape 
similar to the ones considered here can display frozen spa- 
tial chaos [nSp. A numerical search for this manifestation of 
"frozen turbulence" at the level of Navier-Stokes equations is 
currently in progress. We expect this observation to promote 
also experimental work both in the area of pattern formation 
and in hydrodynamics. 



Finally, as an illustrative measure of the asymptotic stabil- 
ity of this solution. Fig. od) displays the value of the Lya- 
punov potential V[ijj] evaluated along the time evolution of the 
field tlj{x,y,t) — A{x , t) sm{Try) . The functional decreases 
in time, confirming t he c onsistency of the STMA with the ex- 
act dynamics of Eq. (2.2), and the potential asymptotically ap- 



proaches a constant value indicating that the field has reached 
a local minimum of V. 
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V. CONCLUSIONS 

We have given evidence of the existence of stationary, sta- 
ble, longitudinally chaotic spatial configurations induced by 
undulated boundaries in a simple two-dimensional reaction- 
diffusion model that does not otherwise display any kind of 
chaos. We have demonstrated that these type of boundaries 
can be convincingly mapped into spatially-peiiodic paramet- 
ric modulations in a one-dimensional approximation to the 
original system. In a dynamical systems approach to the study 
of stationary solutions, these modulations play the role of a 
temporal time-periodic forcing capable to drive a nonlinear 
second-order ODE into chaotic behavior The diffusive char- 
acter of our original model ensures precisely that the relevant 
ODE is in fact, second order and that the presence of chaotic 
stationary solutions is expected at this level of approximation. 
Previous results exist [n,n8Ul^] showing that spatial periodic 
modulation of some parameters intrinsic to the dynamics may 
originate spatial chaos in relatively simple reaction diffusion 
models. However, to our knowledge, the present is the first 
example in which an straightforward dynamical systems ap- 
proach is used to establish the existence of disorder in two- 
dimensional systems due to the influence of the boundaries. 

The consequence of our analysis is that chaotic configu- 
rations should exist in virtually any of the experimental sys- 
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